Nuclear correlations and the r-process 
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We show that long-range correlations for nuclear masses have a significant effect on the synthesis 
of heavy elements by the r-process. As calculated by Delaroche et al. [H, these correlations suppress 
magic number effects associated with minor shells. This impacts the calculated abundances before 
the third r-process peak (at mass number A w 195), where the abundances are low and form 
a trough. This trough and the position of the third abundance peak are strongly affected by 
the masses of nuclei in the transition region between deformed and spherical. Based on different 
astrophysical environments, our results demonstrate that a microscopic theory of nuclear masses 
including correlations naturally smoothens the separation energies, thus reducing the trough and 
improving the agreement with observed solar system abundances. 



The rapid neutron capture process (r-process) leads to 
half of the heavy elements beyond iron. In this nucleosyn- 
thesis process many neutrons are captured by seed nuclei 
(i.e., iron group and nuclei up to A « 90) in time scales 
shorter than beta decay. The pathway for the r-process 
thus involves extremely neutron-rich nuclei far from sta- 
bility and close to the neutron drip-line. However, only 
few of these very exotic nuclei can be produced in current 
rare isotope beam facilities. Even with next-generation 
facilities not all nuclei relevant for the r-process will be 
reached. In view of the experimental difficulties, calcu- 
lations based on theoretical models of the reaction chain 
are essential to understand nucleosynthesis and the origin 
of the heaviest nuclei in the universe. 

Even more than the nuclear physics uncertainties, the 
astrophysical environment of the r-process is far from un- 
derstood [2]. The rapidity and the high neutron densi- 
ties required for this process point to the most neutron- 
rich objects in the universe: neutron stars. Many stud- 
ies have been performed for the high-entropy, neutrino- 
driven wind subsequent to core-collapse supernovae (see 
Q for a review). The conditions necessary for the r- 
process have been well established but are not found in 
recent hydrodynamic simulations [3-6]. The merger of 
two neutron stars was proposed in [7[ and indeed simu- 
lations can reach the extreme conditions required for the 
r-process (see, e.g., (HQ). However, neutron star mergers 
are not occurring early enough in our galaxy to explain 
the heavy r-process elements in very old stars [lOj . 

Typically, the r-process is modeled taking as theoreti- 
cal input some scenario of the astrophysical environment 
together with specific models for nuclei where experimen- 
tal information is not available. While individual models 
can be judged by the quality of agreement for the pre- 
dicted final abundances, the dependence on the nuclear 
physics entering the different models and which physics 
features cause discrepancies is not fully understood. In 



this Letter, we show that some features of the final abun- 
dances can be traced to the treatment of correlations 
in the calculated masses. It has been found that long- 
range correlations associated with quadrupole degrees 
of freedom can significantly improve the description of 
nuclear masses with respect to Hartree-Fock-Bogoliubov 
(HFB) theory, using either a Skyrme interaction [ll[ or 
the Gogny interaction [l[ . It was also found that the cor- 
relation energy enhances magic number effects at double 
shell closures and weakening them elsewhere [12| , a phe- 
nomenon known as "mutually enhanced magicity" [13| . 
The mass table of Ref. [l| contains both the HFB and the 
correlation contribution permitting us to examine its spe- 
cific effect on nucleosynthesis. In both Refs. jl( and [U] 
the HFB theory was extended by the generator coordi- 
nate method to treat quadrupole correlations, and no pa- 
rameters were assumed beyond those going into Skyrme 
or Gogny energy functionals. 

Since there is uncertainty about the astrophysical en- 
vironment of the r-process, we study the effects of long- 
range nuclear correlations for the different scenarios men- 
tioned above. Specifically, we consider trajectories aris- 
ing from hot or cold, high-entropy neutrino-driven winds 
[3| and from neutron star mergers [14( . In all cases we find 
that correlations reduce the trough in the abundances 
before the third r-process peak (at A & 195). Similar 
reductions have been achieved in nucleosynthesis calcu- 
lations based on mass models with shell quenching [l5| . 
Here we show that this trough depends not only on the 
shell closures, but also on the transition from deformed 
to spherical nuclei. 

The most important input to the reaction network 
from the perspective here is the neutron separation en- 
ergy, which crucially affects the (7, n) photodissociation 
reactions and the (n, 7) capture rates. The separation en- 
ergies obtained from the mass tables of Delaroche, el al. 
[l| are shown in Fig. [U The lefthand and righthand pan- 
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FIG. 1. (Color online) Half the two-neutron separation energy S2n/2 for constant proton number as a function of mass number 
A. The lines represent isotopic chains from Z — 30 to Z — 80. The dots mark the r-process path for the hot wind trajectory 
at freeze-out. The left panel corresponds to the nuclear masses of Delaroche et al. |l| without nuclear correlations. In the right 
panel nuclear correlations are included. 



els correspond respectively to calculated masses in HFB 
theory and in their extended theory including correla- 
tions. Plotted are contours at constant proton number 
of half the two- neutron separation energy 52 n , as a func- 
tion of mass number. Delaroche et al. [1] calculated only 
even-even nuclei, so the needed odd-A separation energies 
will be assigned by interpolation. Also, their correlated 
masses are higher in some nuclei and for those masses we 
use the HFB values. With the theoretical S n values, we 
calculate the neutron-capture cross sections with a sim- 
plified Hauser-Feshbach model [Eq. (39) in Ref. 16]. The 
nuclear binding energy differences also affect beta-decay 
rates, but the sensitivity is not as high [ItJ and we use 
the table of [18| for these transitions. 

We point out two features in the calculated two- 
neutron separation energies (Fig. [1]) that have a major 
effect on the r-process abundances, see also [17j and ref- 
erences therein. First, the abrupt drop of S2n around 
A = 130 and A = 195 corresponds to the magic num- 
bers N = 82 and N = 126, respectively. For these nuclei 
with closed shells, neutron-capture cross sections are very 
small and the photodissociation rate is high for typical 
r-process temperatures. As a result the r-process path 
stops at these nuclei and waits for beta decay. It is here 
that matter accumulates and the abundance peaks form, 
as they are observed in the solar system. The second 
characteristic feature is the smoothness of the topogra- 
phy just before the magic numbers. Without correla- 
tions, the evolution of 62 n with A shows a pronounced 
dip just before the N = 126 magic number, associated 
with the transition from deformed to spherical nuclei. 
This behavior is also present in other mass models [IQI — 

Nuclear correlations strongly affect these two features. 



When correlations are included, the neutron shell clo- 
sures are smoother, similarly to what occurs when intro- 
ducing shell quenching [18]. The less pronounced drop of 
S2n at TV = 82 and N = 126 leads to smaller peaks in 
abundances. The other important effect is the smoothing 
of the dip in S 2n just before N = 126 (A « 180). This 
has a big impact on late neutron captures and on the 
formation and evolution of the trough in the abundances 
before the A = 195 peak. 

The evolution of the astrophysical environment is 
shown in Fig. [2] for the three scenarios we consider. First, 
we take two trajectories of neutrino-driven wind simula- 
tions [3] with the entropy artificially increased to pro- 
duce the third r-process peak (i.e., the same trajecto- 
ries as in [l7[). The solid black line represents a cold 
wind evolution where photo-dissociation is negligible and 
the evolution is characterized by a competition between 
neutron captures and beta decays. The other trajec- 
tory, assuming higher temperatures and densities (dashed 
green line), is closer to a classical r-process [22] in which 
(n, 7) — (7, n) equilibrium is achieved. Moreover, we con- 
sider a neutron star merger trajectory (dotted red line) 
from the simulations of Ref. [14] as applied to the r- 
process (see [8]). Notice the fast drop of the density 
(initially very high) and the small increase of the tem- 
perature (at t « 10 -2 — Is) due to the energy released 
by r-process reactions [23[. 

The resulting abundances are shown in Fig. [3l The 
nucleosynthesis calculations start at high temperature (~ 
1 MeV) with nuclear statistical equilibrium. The initial 
evolution until charged-particle freeze-out is calculated 
with a full network of reactions and including nuclei up to 
Z — 63 [13, [H| . The r-process phase which follows (after 
T ~ 0.25 MeV) is modeled with a reduced network [I?], 
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FIG. 2. (Color online) Evolution of temperature and density 
for different astrophysical environments. The trajectories for 
a neutrino-driven wind with cold (solid black line) and hot 
(dashed green line) r-process are based on the simulations 
of 0| . The trajectory for the neutron star merger 0] is given 
by dotted red line. 



l24j that includes the relevant reactions for the r-process 
(neutron capture, photodissociation, beta decay, alpha 
decay, and fission). The masses from Ref. |l| are used in 
this r-process network. 

The abundances are very different for the three trajec- 
tories (see Fig. [3j) because they depend very much on the 
temperature and density evolutions. For example, in the 
hot wind scenario (n, 7) — (7, n) equilibrium is important 
in determining the abundances. In the other scenarios, 
the main competition is between neutron captures and 
beta decays. The red dots in Fig. Q] show the r-process 
path at freeze-out (neutron-to-seed ratio equal one) for 
the hot wind trajectory. The abundances at this point 
thus provide crucial information for the r-process phase. 
Afterwards nuclei beta decay to stability, only changing 
the mass number by beta-delayed neutron emission or 
by capture of the last neutron (see, e.g., [TriE!!)- These 
late-stage reactions lead to the stagger in the resulting 
abundances between peaks. The observed abundances 
are much smoother; the difference may be attributed to 
the approximation used for the neutron captures [ItJ- 
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FIG. 3. (Color online) Abundances with and without nuclear 
correlations based on hot and cold wind and neutron star 
merger trajectories, from the top to the bottom. The black 
dots are solar system r-process abundances. 

None of the scenarios shown here fit the entire range of 
solar system abundances. Very likely, the observed abun- 
dances arise from a superposition of trajectories that in- 
dividually emphasize the lighter or the heavier mass num- 
bers. Our goal here was simply to show relative changes 
due to correlations, in scenarios that lead to the heavier 
elements. 

From Fig. [3] one sees that nuclear correlations affect 
mainly the abundances in the region before the third 
peak (180 < A < 190). Here, the three trajectories show 
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a similar behavior; the correlations reduce the trough. 
However, the depth of the trough and the impact of cor- 
relations depends on the trajectory. For the hot wind 
(top panel), (n, 7) — (7, n) equilibrium is reached and the 
r-process path is given by the neutron separation energy. 
Therefore, the features in S2n become crucial to under- 
stand the final abundances. The abundances at freeze- 
out already show a pronounced trough that is reduced 
and shifted toward higher A during the decay to sta- 
bility. The evolution and fill-up of the trough therefore 
depends on neutron captures. These can be enhanced 
due to the features shown in 52 n /2 ~ 3 — 4 MeV at 
A & 180 without correlations (left panel Fig. [I]). For the 
cold wind (middle panel) the trough and the second peak 
are smaller than for the hot wind. Because of the lower 
temperatures during the r-process in this case the pho- 
todissociation is almost negligible. Therefore, the impact 
of the masses is smaller as it enters mainly through the 
neutron captures. The merger trajectory (bottom panel) 
presents extreme abundances with a very narrow third 
peak. The dip in S2n increases the neutron captures in 
the region before the third peak and moves matter to- 
wards higher mass numbers. This results in a narrower 
trough for the calculation with correlations. 

In summary, correlations smooth out the nucleus-to- 
nucleus variations of nuclear separation energies. This is 
especially pronounced for the region before N = 126 as 
shown by the two neutron separation energy. The dip 
observed in S2n without correlations is also present in 
other mass models Jl9H2l| (see Fig. 6 of [17]). In the 
classical r-process [22) \ the problem of the large trough 
before peaks was previously overcome by quenching the 
shell gap [TBI. The masses of nuclei in the transition 
region (180 < A < 190) between deformed and spheri- 
cal strongly affect the calculated abundances determining 
the shape and position of the third peak and the trough 
before this. We have shown that a microscopic theory 
of nuclear masses including correlations is sufficient to 
smoothen the separation energies and thus reduce the 
trough in the abundances. Our study here shows the 
strong impact of correlations on r-process nucleosynthe- 
sis and motivates further theoretical work to treat the 
nuclear physics as realistically as possible. For the fu- 
ture, we plan to study nucleosynthesis using globally cal- 
culated beta decays and neutron captures consistently 
with the theoretical separation energies. 

We thank F.-K. Thielemann, G. Martmez-Pinedo, and 
T.R. Rodriguez for discussions and comments on the 
manuscript. A. A. is supported by a Feodor Lynen Fellow- 
ship of the Alexander von Humboldt Foundation and by 
the Swiss National Science Foundation. G.F.B. acknowl- 
edges support from the US Department of Energy under 
grant DE-FG02-00ER41132. We also thank the Institute 
for Nuclear Theory, where this work was initiated. 



[2] 
[3] 
[4] 
[5] 

[e; 

[T 

[8 
[9 
[10 

[11 
[12 
[13 
[14 

[15 

[16 

[17 
[18 

[19 
[20; 
[21 
[22; 
[23 

[24 

[25; 



la.arc ones@unib as.chl 
J.-P. Delaroche, M. Girod, J. 
S. Hilaire, S. Peru, N. Pillet 
|Phys. Rev. C 81, 014303 (2010)| 
M. Arnould, S. Goriely , 
|Phys. Repts. 450, 97 ( 2007) 
A. Arcones, H.-T. Janka, 



Libert, 
and G. 



H. Goutte, 
F. Bertsch, 



and K. Takahashi, 



and L. S check, 



Astro n. fe Astrophys. 4 67, 1227 (2007)] 
L. Hud epohl, B. Miiller, H. Janka, A. Marek and G. G. 
Raffelt, |Phys. Rev. Lett. 104, 251101 (2010 )] 

A. Mezza- 
Liebendorfer, 



T. Fischer, S. C. Whitehouse, 
cappa, F. Thielemann, and M. 



Aston. & Astrophys. 517, A80 (2010) 



Vrcones 



and 



H.-T. 



Janka, 

|Astron. fe Astrophys. 526, A160 (2011)| 

J. M. Lattimer and D. N. Schramm, 

|Astrophys. J. 192, L145 (1974) | 

C. Freiburghaus, S. Rosswog, and F.-K. Thielemann, 
|Astrophys. J. 525, L121 (1999) | 

S. Goriely, A. Bauswein, and H.-T. Janka, 
|Astrophys. J. 738, L32 (201 1)| 

D. Ar gast, M. Samland, F. -K. Thielemann, and Y.-Z. 
Qian, Astron. fe Astrophys. 416, 997 (2004) | 

M. Bender, G. F. Bertsch, and P.-H. Heenen, 
|Phys. Rev. C 73, 034322 (2006) | 

M. Bender, G. F. Bertsch, and P.-H. Heenen, 
|Phys. Rev. C 78, 054312 (2008) | 

D. Lunney, J. M. Pearson, and C. Thibault, 
|Rev. Mod. Phys. 75, 1021 (2003) | 

S. Rosswog, M. Liebendorfer, F.-K. Thielemann, M. B. 
Davies, W. Benz, and T. Piran, Astron. & Astrophys. 
341, 499 (1999). 

B. Chen, J. Dobaczewski, K.-L. Kratz, K. Lan- 
ganke, B. Pfeiffer, F.-K. Thielemann, and P. Vogel, 
|Phys. Lett. B 355, 37 (1995)| 

S. E. Woosley, W. A. Fowler, J. A. Holmes, and B. A. 
Zimmerman, Astrophysical Reaction Rates, Preprint 
OAP-422 (California Institute of Technology, W. K. Kel- 
logg Radiation Laboratory, 1975). 

A. Arcones and G. Martmez-Pinedo, 
|Phys. Rev. C 83, 04 5809 (2011)1 

P. Moller, B. Pfeiffer, and K.-L. Kratz, 

Phys. Rev. C 67, 055802 (2003), see Table a t 

http://tl6web.lanl.gov/Moller/publications/tpnff.data 

P. Moller, J. R. Nix, W. D. Myers, and W. J. Swiatecki, 

|At. Data Nucl. Data Tables 59, 185 (1995)| 

J. M. Pearson, R. C. Nayak, and S. Goriely, 

|Phys. Lett. B 387, 455 (1996)] 

S. Goriely, N. Chamel, and J. M. Pearson, 
|Phys. Rev. C 102, 152503 (2009) | 

K. Kratz, J. Bitouzet, F. Thielemann, P. Moeller, and 

B. Pfeiffer, Astrophys. J. 403, 216 (1993). 

B. D. Metzger, A. Arcones, 

E. Quataert, and G. Martmez-Pinedo, 
|Monthly Notices of the Royal Astronomical Society 402, 2771 (2010 

C. Freiburghaus, J.-F. Rembges, T. Rauscher, E. Kolbe, 

F. -K. Thielemann, K.-L. Kratz, B. Pfeiffer, and J. J. 
Cowan, Astrophys. J. 516, 381 (1999). 

M. Mumpower, G. McLaughlin, and R. Surman, 
larXiv:1109.3"613l 



